Logarithmic temperature profiles in turbulent Rayleigh-Benard convection 
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We report resuhs for the temperature profiles of turbulent Rayleigh-Benard convection (RBC) in 
the interior of a cylindrical sample of aspect ratio T = D/L — 0.50 {D and L are the diameter and 
height respectively). Results from experiment over the Rayleigh number range 4 x 10^^ ^ Ra < lO'^^ 
for a Prandtl number Pr ~ 0.8 and from direct numerical simulation (DNS) at Ra = 2 x 10^^ for 
Pr = 0.7 are presented. We find that the temperature varies as A * ln{z/L) + B where z is the 
distance from the bottom or top plate. This is the case in the classical as well as in the ultimate 
state of RBC. From DNS we find that A in the classical state decreases in the radial direction as 
the distance from the side wall increases and becomes small near the sample center. 

PACS numbers: 47.27.te,47.32.Ef,47.20.Bp,47.27.ek 



Turbulent convection of a fluid contained between two 
horizontal plates separated by a distance L and heated 
from below (Rayleigh-Benard convection or RBC) [iHSj is 
a system in fluid mechanics with many features that are 
of fundamental interest. It is also a phenomenon with nu- 
merous astrophysical geophysical 0413, ^'^'^ tech- 
nological 13] applications. Nonetheless some of its prop- 
erties remain poorly explored and understood. A "clas- 
sical" state of RBC exists below a transition range to an 
"ultimate" state; the transition range extends over more 



than a decade from Ra]' to Ra2 (Ra is a dimension- 
less measure of the applied temperature difference). For 
simplicity we shall characterize this range by Ra* which, 
for the parameters of our work, is about 10^'' [3 [l5| . 
For the classical state it is known from experiment (see, 
for instance, 16-23) that approximately half of the ap- 
plied temperature difference AT is sustained by two thin 
thermal boundary layers (BLs), one just below the top 
and the other just above the bottom plate. These BLs 
are of Prandtl-Blasius type, i.e. laminar, albeit fluctu- 
ating 2J, |25|. The entire interior of the sample, known 



as the "bulk" , is then approximately isothermal in the 
time average, but it also undergoes vigorous local tem- 
perature fluctuations. At a more detailed level it has 
long been recognized that the bulk actually sustains small 
temperature gradients, but these gradients were believed 
to be more or less independent of vertical position and 
the total temperature drop across the bulk is known to be 
much smaller than that across the BLs (see, for instance. 



in the experiment, the temperature T{z) and its root- 
mean-square (rms) fluctuations cr{z) vary logarithmically 
as a function of the distance z from the bottom plate. 
While a logarithmic dependence was predicted for the ul- 
timate state above Ra* where the BLs are turbulent and 
believed to extend throughout the entire sample [i^, to 
our knowledge there is no theory at present that predicts 
the bulk properties for Ra < Ra*. We believe that the 
discovery of logarithmic profiles is an important step to- 
ward developing a more fundamental understanding of 
the bulk. On the one hand, its origin may be sought in 
the diffusion of the excess or deficit of enthalpy carried 
by plumes. On the other hand, the logarithmic depen- 
dence suggests a relationship to the well known logarith- 
mic variation of the velocity in turbulent shear flows dis- 
cussed originally by von Karman and Prandtl [13 
(for a recent review, see [3l|). 

We present experimental measurements for a Prandtl 
number Pr ~ 0.8 over the range 4 x 10^^ ^ Ra < 10^^ and 
analyze direct numerical simulations (DNS) for Pr = 0.7 
and Ra = 2 x 10^^ [13 for a cylindrical sample of 
aspect ratio F = D/L = 0.50 {D is the diameter). 
Both experiment and DNS show that, through most of 
the bulk, the dimensionless time-averaged temperature 
e(z) = [(T(z)) - T^]/iTb - Tt) (we denote the time av- 
erage by (...) and = (Tb + Tt)/2 with Tfc and Tt the 
temperatures at the top and bottom of the sample) can 
be represented well by 



e(z) ^ A*ln{z/L) + B 



(1) 



We find that, beyond a thin boundary layer unresolved 



We also find that the root-mean-square temperature flue- 
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tuations cr(z) = {[T{z)- {T{z))fY/'^ / {Tb-Tt) are consis- 
tent with a logarithmic dependence on z, and represent 
them by 



a{z) ^C*ln{z/L) + D 



(2) 



The DNS data show that the amphtude A{r) is largest 
near the side wall and decreases gradually as the distance 
R — r from the wall increases (i? = D/2 is the sample 
radius and r the radial coordinate) . At the sample center 
A nearly vanishes. The experimental data cover a wide 
range of Ra but are for a single {R — r)/L = 0.0045. This 
location is, however, well inside the bulk; the DNS data 
at Ra = 2 X 10^^ showed that the viscous BL only extends 
to (i?-r)/Lw 0.0008. 

The measurements were made with a large cylindrical 
sample of height L = 2.24 m and diameter D = 1.12 
m known as the High-Pressure Convection Facility II 
(HPCF-II) which was placed in an even larger pressure 
vessel known as the "Uboot of Gottingen" at the Max 
Planck Institute for Dynamics and Self Organization in 
Gottingen, Germany [33[. The Uboot and HPCF-II were 
filled with the gas sulfur hexafluoride (SFg). During the 
measurements the HPCF-II was completely sealed. The 
Prandtl number Pr = v / k. {u is the kinematic viscos- 
ity and n the thermal diffusivity) was 0.79 (0.86) near 
Ra = 4 X 10^^ (lO'^^). The measurements were made at 
Tm — 21°C and at various pressures up to 19 bars. The 
Rayleigh number is given by Ra = agt^Tl? j kv. Here 
the isobaric thermal expansion coefficient a, as well as k 
and were evaluated at T^, and g is the acceleration of 
gravity. Typically AT had values in the range from 4 to 
16 K. 

The sample was tilted slightly, with its axis at an an- 
gle of 14 mrad relative to gravity. We ensured that the 
tilting had no effect on our results. This is understand- 
able because at these high Ra values no pronounced large 
scale circulation (LSC) exists IJ]. Two sets of ther- 
mistors were installed for the temperature-profile mea- 
surements. One was located at what would be the pre- 
ferred down-flow orientation at lower i?a, and the other 
was removed from the first in the azimuthal direction 
by an angle tt. Each set consisted of eight thermome- 
ters which were located in the fluid 1.0 ± 0.1 cm from 
the side wall, i.e.. at a radial position [R — r)/L — 
0.00445 ± 0.0005. The eight thermistors were located at 
z = 4.0, 6.1, 8.1, 12.1, 16.1, 32.2, 64.2, and 110.5 cm, with 
an uncertainty of the vertical position of 0.1 cm. 

Two typical measurements of the sixteen time- 
averaged temperatures are shown in Fig. [1^ as a function 
of z/L on a logarithmic scale. The top (bottom) one is 
for Ra = 1.18 x lO^^ g4 ^ iQi^) in the classical (ulti- 
mate) state. For each case one sees that the data sets at 
the two azimuthal orientations differ slightly from each 
other. We attribute this to the influence of remnants of 
the LSC on the temperature profiles f33|. To compen- 
sate for this effect, we henceforth consider the average at 
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FIG. 1. The measurements are for a radial position [R — 
r)/L = 0.0045. (a): The time-averaged temperature Q{z) = 
[{T{z)) — Tm]/{Tb — Tt) as a function of the vertical position 
z/L. The vertical dotted line marks the sample center at 
zjL = 0.5. The open squares (diamonds) are results at the 
preferred down-flow (up-flow) orientation of the LSC. The 
solid symbols are their averages. The solid squares (red on- 
line, z/ L 0.08) were used for the fits of Eq. [T]to the data. 
The solid lines are those fits. The upper data set is in the 
classical regime for Ra = 1.18 x 10^^ and Pr — 0.787. It gave 
A = -0.0162 ± 0.0004 and B = -0.00218 ± 0.0007 where the 
uncertainties are the standard errors obtained from the fit. 
The dashed line is a fit of a power law Q{z) — A* (zjV)'' -\- B 
to the data with z/L < 0.08 which yielded C = -0.09 ± 0.05. 
The lower data set is for the ultimate state at Ra = 9.64 x 10^* 
and Pr = 0.862. It gave A = -0.0082 ± 0.00016 and 
B = -0.0352 ± 0.0005. (b): The root-mean-square fluc- 
tuations a{z) = {[T{z) - {T{z))fy^ /{n - Tt) as a func- 
tion of z/L corresponding to the data in (a). The sym- 
bols are as in (a), and the fit parameters of Eq. [2] are 
C = -0.00192±0.00013 and D = 0.00565±0.00044 for the up- 
per (classical) set at Ra = 1.18 x 10^^ and C = 0.0019±0.0002 
and D = 0.0030 ± 0.0007 for the lower (ultimate) set at 
Ra = 9.64 X 10^*. 



each vertical position of the two data sets, as shown by 
the solid symbols in the figure. Except at the largest z/L, 
these data fall on straight lines and thus are represented 
well by Eq. [T] for more than a decade. 

In a sample that conforms perfectly to the Boussi- 
nesq approximation we would expect another logarith- 
mic dependence emanating from the top plate to meet 
the data shown in the figure at Tm {i.e. at O = 0) 
and z/L = 1/2. However, in the experiment we find 
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that e(z/L = 1/2) < 0, albeit only by 0.006 (0.028) for 
Ra= 1.18x10" (9.64x10^^). We do not know the reason 
for this offset. However, it will necessarily lead to a small 
departure from the logarithmic dependence because the 
two branches, one coming from the bottom and the other 
from the top plate, must form an analytic function with 
a continuous derivative at z/L =1/2 where they meet. 
For a quantitative analysis we therefore fit Eq. [T] only to 
the five points with z/L 0.08. The resulting functions 
are shown as the solid lines in the figure. 

It should be mentioned that, as for any log-behavior 
in the limited range of 1.5 decades, one can also ob- 
tain a good fit to the data with the power law 8(z) = 
A * (z/L)'' + B. This is shown by the dashed line in 
Fig.[T^. However, the resulting exponent C, = 0.09 ±0.05 
is quite small and unlike any other known exponent rel- 
evant to this system 0]. Further, its uncertainty is 
not much smaller than its value. Thus we shall continue 
the presentation of the results in terms of the logarithmic 
function Eq. [TJ 

The rms temperature fluctuations a are shown in 
Fig. [T}d. They too are described well by the logarith- 
mic form. Also in this case the relevant equation (Eq. [5]) 
was fit to the data only for z/L < 0.08 to determine C 
and D. 

The parameters A, B, C, and D, as well as the devi- 
ation of the center temperature Tc = T{z/L = 0.50) 
from Tm, are shown in Fig. [5] as a function of Ra. As 
had been reported before for measurements of the Nus- 
selt and the Reynolds number [IJ] and reported above. 



there is a range of Ra which extends from Raj to Ra2 
(the vertical dotted lines in the figure) over which the 
transition from the classical to the ultimate state takes 
place. The locations of Raj and Ra2 are particularly no- 
ticeable in the data for B and $. In the transition region 
the parameters scatter much more than above or below it 
because the state assumed by the system can vary from 
one experimental point to another. 

In Fig. [3] we show the results for 6(^) obtained from 
DNS at Ra = 2 X lO^^ ^nd Pr = 0.7^. Those profiles 
are for {R — r)/L — 0.0045, i.e. the same radial position 
as that of the experiment, and are based on azimuthally 
and time averaged temperature data. In this (and the fol- 
lowing) figure we show the measured profiles in the sim- 
ulations both in the top and the bottom half of the sam- 
ple. The difference between these two is an indication of 
the statistical uncertainty due to limited time averaging, 
which of course is a more serious limitation in numerics 
than it is in the experiment. Just as for the experimental 
data, the profile at larger z/L can be described well by 
Eq. (m . Figure |4] gives the DNS results for A and B as 
a function of the radial position {R — r)/L, based on the 
temperature data for 10~^ < z/L < 10^^. In this figure 
one sees that there is excellent agreement between the 
values of A and B (when corrected for the temperature 
measured a.t z/L = 0.50 (Fig. ^)) measured in the ex- 
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FIG. 2. The parameters A, B, C, and D obtained by fitting 
Eqs. [1] and [2] to the experimental temperature and fluctua- 
tion profiles, and the deviation </> = (Tc — Tm)/AT of the 
temperature Tc = T{z/L = 1/2) from the mean temperature 
Tm, all as a function of Ra. All data are for a radial position 
{R — r)/L — 0.0045. The vertical dotted lines indicate the 
locations of Ra^ and Ra2. The solid lines (red online) are 
fits of the function A = a * logio(Ra.) + & to the data with 
Ra < Raj. The extrapolation to Ra = 2 x 10^^ (open cir- 
cles, green online) yielded A = -0.0212, B = -0.0296, and 
<E> = -0.0098. 



periment and simulations. In addition, the figure reveals 
that the magnitude of A is largest near the side wall and 
that it decreases (approximately logarithmically) as the 
sample interior is approached. Although it remains fi- 
nite even at the sample centerline [{R — r)/L — 1/4], one 
might conjecture that it would indeed vanish in a sample 
of larger aspect ratio. Thus the DNS indicates that the 
logarithmic vertical temperature profile is related to the 
existence of the side wall; however, it penetrates radially 
deep into the sample, over a distance that is well over an 
order of magnitude larger than the BL thickness. 
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stant or varying linearly and slowly in space, we find that 
Q varies logarithmically with distance from the plates 
over a wide range of z/L. The root-mean-square tem- 
perature fluctuations show similar variations. The ampli- 
tude of the logarithmic profile is largest near the side wall 
and becomes small as the sample center is approached. 
The origin of the logarithmic profile remains unclear. On 
the one hand one may speculate that it is the result of the 
diffusion of enthalpy carried from the BLs into the inte- 
rior by plumes; but a model for this process which would 
yield a logarithmic distribution is not known to us. On 
the other hand, the logarithmic variation suggests a pos- 
sible relationship to the well known logarithmic velocity 
profiles in shear flows 29l-31|. 

In the ultimate state, which exists near Ra = lO^'^ 
range lO"'' < z/L < 0.1 to the average of the two profiles (the and above [H [H, it was predicted [28] that the BLs 
fitting interval is indicated by the two short dashed vertical 
lines). 



FIG. 3. The temperature 0{z) as a function of z/L at a 
radial position {R -r)/L = 0.0045 horn DNS for Ra = 2 x 
10^2, Pr = 0.7, and V = 1/2. The solid red (dashed blue) line 
indicates the profile measured from the bottom (top) plate as 
a function of z/L (here z is taken as the distance from the 
nearest plate). The thin line represents a fit of Eq.[T]over the 
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FIG. 4. DNS results for A and B in Eq. [T] as a function of 
{R~r)/L on a logarithmic scale. The solid red (dashed blue) 
lines indicate the result obtained above the bottom (below 
the top) plate. The open circles (green online) correspond to 
the extrapolations of the experimental data for A and B as 
shown in Figs. [2^ and b to Ra = 2 x 10^^. The open diamond 
indicates the value of B corrected for $, see Fig. [2^. 



are turbulent and that they extend vertically throughout 
the entire sample; thus there is no "bulk" in the same 
sense as there is for the classical state. A logarithmic 
temperature profile due to the turbulent BLs was pre- 
dicted to extend from each plate deep into the sample, 
with the two profiles meeting at half height. Indeed, the 
experimental measurements in the ultimate state do find 
a logarithmic dependence of the temperature on the ver- 
tical coordinate. Unfortunately, these large values of Ra 
are not yet accessible to DNS (and will not be for some 
time), and experimental results are available only for one 
radial position. Thus it is not known whether the loga- 
rithmic variation persist throughout the sample, as one 
would expect on the basis of the prediction. 
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In this Letter we reported on results obtained by using 
a combination of experiment and DNS to study the in- 
terior of turbulent RBC. For the classical state, which 
exists for Ra < Ra* ~ 10^^ and which has laminar 
BLs adjacent to the top and bottom plates, we find 
that the bulk which is found between these two layers 
sustains a non-trivial and interesting temperature field 
Q{z/ L,r / L). Whereas it had generally been assumed 
that the temperature in the sample interior is either con- 
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